Quantum phase transitions of light. 
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Recently, condensed matter and atomic exper- 
iments have reached a length-scale and tempera- 
ture regime where new quantum collective phe- 
nomena emerge. Finding such physics in systems 
of photons, however, is problematic, as photons 
typically do not interact with each other and can 
be created or destroyed at will. Here, we intro- 
duce a physical system of photons that exhibits 
strongly correlated dynamics on a meso-scale. By 
adding photons to a two-dimensional array of 
coupled optical cavities each containing a single 
two-level atom in the photon-blockade regime, 
we form dressed states, or polaritons, that are 
both long-lived and strongly interacting. Our 
zero temperature results predict that this pho- 
tonic system will undergo a characteristic Mott 
insulator (excitations localised on each site) to 
superfluid (excitations delocalised across the lat- 
tice) quantum phase transition. Each cavity's im- 
pressive photon out-coupling potential may lead 
to actual devices based on these quantum many- 
body effects, as well as observable, tunable quan- 
tum simulators. 

The Jaynes-Cummings model is arguably the most 
important model for understanding light-matter inter- 
actions. It describes the interaction of a single, quasi- 
resonant optical cavity field with a two-level atom. The 
coupling between the atom and the photons leads to op- 
tical nonlinearities and an effective photon-photon repul- 
sion. Perhaps the most extreme demonstration of this 
photonic repulsion is photon blockade, demonstrated re- 
cently by Birnbaum et al. 0, where photonic repulsion 
prevents more than one photon from being in the cavity 
at any one time. Photon blockade was initially theoreti- 
cally described with a four-state system ^ , with multipli- 
cation of the weak Kerr nonlinearity effected by placing 
a large number of atoms within each cavity. However, 
it was quickly realised that the photonic blockade mech- 
anism does not persist in the limit of many atoms 0, 
rapidly degrading as the number of atoms per cavity is 
increased Later Rebic et al. showed that the non- 
linear interaction afforded by placing a single two-level 
atom inside a cavity would suffice for realising photon 
blockade jS]. This observation was highly significant as 
it allowed the full weight of the Jaynes-Cummings model 
to be used to attack and understand this problem. 



To create an atom-photon system whose dynamics mir- 
ror those traditionally associated with strongly inter- 
acting condensed matter systems, we consider a two- 
dimensional array of photonic bandgap cavities. Each 
cavity contains a single two-level atom, quasi-resonant 
with the cavity mode. Evanescent coupling between the 
cavities due to their proximity allows inter-cavity photon 
hopping. This configuration is depicted schematically in 
Fig. n^a), where we have explicitly chosen three near- 
est neighbours per cavity (coordination number z — 3), 
for reasons explained below. Because we are considering 
small cavities, with volumes of order where A is the 
wavelength of the light, there will be strong atom-photon 
couplings that will dominate over the spontaneous emis- 
sion. Fig. n^b) shows the pertinent energy scales within 
one of the cavities. In the classical limit without two-level 
atoms, an array of coupled photonic bandgap cavities has 
been described for novel waveguiding applications 0, H| , 
and in the quantum regime a two cavity arrangement has 
been proposed as a Q-switch .Q] • We consider here for the 
first time, the rich dynamics of a two dimensional lattice 
of quantum cavities. 

Interest in interacting boson systems grew out of work 
on the metal-insulator transition in Fermi systems |inl |. 
The seminal paper of Fisher et al. employed an infi- 
nite range mean-field theory to qualitatively describe the 
superfluid-insulator phase transitions of a bosonic gas in 
random and periodic potentials, with and without disor- 
der. Recently, it has been shown that this Bose-Hubbard 
model can be realized in cold atom optical lattices p^ . 
The Bose-Hubbard Hamiltonian seeks to describe the 
many-body dynamics of strongly interacting bosons, 
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where i,j range over all sites in the lattice (fixed in our 
case) , tij is the hopping energy of bosons between sites i 
and J, U is the on-site repulsion between particles, and 
fli is the bosonic annihilation operator on site i. In our 
photonic-atom superlattice, the Jaynes-Cummings inter- 
action provides the effective on-site repulsion, and the 
hopping term is represented by the evanescent coupling 
between the cavities. However, our physical system is 
fundamentally different from the Bose-Hubbard case: in 
particular, the effective U is not constant in the Jaynes- 
Cummings system, but decreases as the number of pho- 
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tons in each cavity goes up; and the conserved particles 
here are not pure bosons, but dressed photons or polari- 
tons, which are a mixture of the spinor atom and pho- 
tons in each cavity. Thus, the hopping term corresponds 
to the overlap of these extended composite photon-atom 
states, which are the conserved particles in the model, 
not the photons. 

To motivate the search for Hubbard model type inter- 
actions within a superlattice of photonic bandgap cavi- 
ties, we first discuss the Hamiltonian for a single two-level 
atom in a quasi-resonant cavity, 

H''^ ^ ea+a^ +uja'^a + P {a+a + a-a'') , (2) 

where cr+ and (T_ (a, a^) correspond to the atomic (pho- 
tonic) raising and lowering operators, respectively. The 
transition energy of the atomic system is e, the cavity 
resonance is uj and the cavity mediated atom-photon cou- 
pling is P, which is implicitly assumed to be real for our 
purposes. Defining the atomic states as \g) and |e) for 
ground and excited state respectively, we introduce the 
detuning, A = lu — e. The eigenstates of Eq. |21are the 
dressed states which we label as |±,n) defined in 
the Methods section, where n is the number of excita- 
tions in the cavity. The ground state for the Jaynes- 
Cummings Hamiltonian is qualitatively different from the 
other dressed states, implying a non-trivial form of the 
raising operator 14']. This constitutes a further depar- 
ture from usual Hubbard-like condensed-matter models, 
where the raising operator is not dependent on the num- 
ber of excitations. 

In Fig.[21we show some of the eigenvalues of the Jaynes- 
Cummings system with the photon energy subtracted 
for ease of comparison, i.e. we are plotting i?|±^ra) ~ i^^- 
The eigenenergies are well known and are i?|±,„) = 

nu! ± x{^) — A/2 where x("-) = V^^/3^ + A^/4 is the n 
photon generalised Rabi frequency. Displaying the eigen- 
values in this way immediately allows us to connect the 
Jaynes-Cummings Hamiltonian with the on-site repul- 
sion in the Hubbard Hamiltonian. The on-site repulsion 
is evinced by the increasing energy separation with n. 

The Hamiltonian for our extended Hubbard-like sys- 
tem is given by a combination of the Jaynes-Cummings 
Hamiltonian with photon hopping between cavities and 
the chemical potential term, 

H = ^ n^^' + («Iaj + aja,;) - ^ fiM, (3) 

where the inter-cavity hopping occurs with frequency 
for nearest neighbours, and n ~Q otherwise, Ni is 
the total number of atomic and photonic excitations (the 
conserved quantity in our system) and /i^ is the chemical 
potential at site i in the grand canonical ensemble. For 
our proof of concept calculation we assume zero disorder 
and = (1 for all sites. To most effectively explore the 
important regime where on-site repulsion dominates over 



hopping, we consider a lattice with as few nearest neigh- 
bours as possible to achieve a two-dimensional network, 
that is three nearest neighbours. Altering the number 
of nearest neighbours does not qualitatively affect our 
results. 

To gain insight over the properties of the full Hamilto- 
nian, we will employ a mean-field approximation. Mean- 
field theories [1^ lla ll^ Si'^s good qualitative and quan- 
titative descriptions of these systems, comparing well to 
Monte Carlo simulations [lil fl^ . These mean- field ap- 
proaches have also been extended to dipolar bosons [23 
and boson-fermion atomic mixtures 0, , as well as 
to the theory of exciton and exciton-polariton (electron- 
hole pair plus photon) condensation We introduce 
a super-fluid order parameter, = (ui), and employ 
the decoupling approximation ajaj — {a\)aj -f {aj)a\ — 
{a\){aj). The resulting mean- field Hamiltonian can be 
written as a sum over single sites, 

i 

where z = 3 is the number of nearest neighbours. To 
obtain the system's zero temperature properties, we use 
the procedure of Refs. llTl which is outlined in the 
Methods section. When -0 = we have a Mott phase, 
characterised by a fixed number of excitations per site 
with no fluctuations, and "0^0 indicates a superfluid 
phase. The boundary between the '0 = and 
phases denotes where a quantum phase transition in this 
system will occur: a quantum phase transition of light. 
In general, we expect that when photon-photon coupling 
(on-site repulsion) dominates over hopping, the system 
should be in a Mott phase, and when the converse is 
true, the system will be in a superfluid phase. 

We can determine the extent of the Mott lobes in the 
limit of very small hopping from inspection of Eqs.[5]and 
0] First note that with our definition, £'|-,n) < ^\+,n) i w6 
only need consider the negative branch for the purposes 
of determining the ground state. Furthermore, a change 
in the total number of excitations per site will occur when 
E\-.n+i) ^ fJ-in + 1) = E\_ n) ^ We can determine 
the critical chemical potential, iidn), where the system 
will change from n to n -I- 1 excitations per site as 

Hc{n) ^ UJ ~ [x{n + I) - xin)] ■ (5) 

These boundary chemical potentials are shown in Fig. |31 
as a function of detuning in the limit of very small k, and 
these also account for the asymptotes in Fig.0] In Fig.^ 
we show the complete phase diagram of the mean field 
solution. Because of the large parameter space, there 
are several cases that we consider. Note that our results 
are only shown for /i < a; as in the limit ^ > lu the 
minimisation method does not converge [/Xc(oo) = w]. 
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although in this Umit, the quantum ground state should 
correspond to a coherent state of excitations (superfluid 
phase). 

The dynamics illustrated in Fig. 0] are extremely rich. 
Hubbard-like dynamics can be seen in the three sub- 
figures of Fig. ^ which show t/; as a function of k and fi for 
three different detunings. In this plot, ip = corresponds 
to stable Mott lobes, with the number of excitations in- 
creasing with fi. The regions to the right correspond to 
/i 7^ 0, and in these regions the system will be found in 
a superfluid phase, i.e. the stable ground state at each 
site corresponds to a coherent state of excitations over 
the \ — ,n) branch. The size of the Mott lobes varies with 
A, with the largest Mott lobes found on resonance. 

To confirm the number of photons in each Mott lobe, 
we show in Fig. [3 the average number of excitations per 
site in the grand-canonical ensemble, p = ~dEg(ip — 
)/dfi, as a function of k and ^ for A = 0. These 
plataues indicate regions of constant density and incom- 
presibility, both characteristic of the gapped Mott insu- 
lator phase. 

POTENTIAL IMPLEMENTATIONS 

Until now, we have deliberately avoided specific im- 
plementations, as the above analysis should be applica- 
ble to all atom-photon systems with strong coupling. In 
this section we specifically consider some candidate im- 
plementations that may be realisable with present-day or 
near-term technology. 

One candidate system would be a photonic bandgap 
structure in diamond, with individual negatively-charged 
nitrogen- vacancy (NV) centres placed in each cavity. Mi- 
cromachining of diamond is presently being explored for 
similar purposes, see 0,113 1 and modelling has shown 
that suitable cavities are in principle fabricatable [j^ . 
The application of diamond for quantum computing and 
quantum optics applications has recently been reviewed 
|27| . For NV diamond, the two-state transition of inter- 
est is at 637 nm ~ 3 x 10^^ Hz. Calculations presented 
m Ref.H showed that A single-mode photonic bandgap 
cavities would have /3 ~ lO^'^Hz. Assuming that photon 
hopping limits the cavity Q, we can approximate tunnel- 
ing frequency by k = ^IQ, so for n < (3 would require 
Q > 10^. This is demanding, but we note that Q ~ 10^ 
has been achieved in silicon-on-silica photonic-bandgap 
cavities [l^, although we stress that complete modeling 
is necessary to precisely determine the required geom- 
etry. If the diamond substrate is ultra-purity Type Ila 
(low nitrogen) diamond, individual NV centres can be im- 
planted using single-ion implantation techniques |27Ll29l| . 
Finally we note that the resonance frequencies for the 
photonic bandgap cavities will be extremely difficult to 
tune post-creation, however the Stark effect can be used 
to tune the NV centres as required 0| to allow an ex- 



ploration (either statically or dynamically) of the phase 
space shown in Fig. 0] 

Another candidate system where such effects may be 
observable is the microwave strip line resonator for cir- 
cuit QED recently demonstrated by Walraff et al. [3lj |. 
In such systems, the effective mode volume of the cav- 
ity can be much less than A*^, with 10~^A'^ demonstrated 
in Ref. |31J, with the rate of loss of photons out of the 
cavity (analogous to the k described in this paper) being 
~ 10'^ Hz. When such values are combined with their 
atom-photon couplings, which were of order (3 - lO'^ Hz, 
then the realised cavities would seem to be well-suited 
to observing the effects we have described. One advan- 
tage of the strip lines over the photonic bandgap case, 
however, would be that generating structures with con- 
trolled (and perhaps varying) numbers of nearest neigh- 
bours would be much easier. 

The above two examples do not constitute the limit of 
possible implementations. It may be possible to achieve 
these phase transitions using Rydberg atoms in high- 
Q superconducting cavities. Such systems have shown 
impressive results including generation of one and two- 
photon states in the cavity 1331 ■ An alternative archi- 
tecture would be an array of microcavities containing 
single atoms, with the cavities connected by optical fi- 
bres 01 . Again, this system would appear to offer the 
desirable properties of high atom-photon coupling, with 
relatively low hopping terms, and classical routing would 
permit interesting lattice structures. Optical quantum 
dots, for example In As structural quantum dots in a pho- 
tonic bandgap lattice, may also afford interesting possi- 
bilities for the realisation of these (or similar) quantum 
phase transitions. 

CONCLUSIONS 

We have shown how a Mott insulator to superfluid 
transition can be realized with atom-dressed photons in 
a cavity superlattice. Calculations for zero temperature 
give a phase diagram analogous to, but fundamentally 
different from, the Bose-Hubbard model. 

If one could Q-switch the cavities, for example by 
bringing each cavity into close proximity to a near field 
probe, it may be possible to out-couple all of the cavities 
simultaneously. This would generate a self-ordered, two- 
dimensional array of single photon emitters, which could 
be extremely beneficial for quantum information applica- 
tions. Because there is significant flexibility in construct- 
ing exotic lattice geometries, one may be able to engineer 
an extended Bose-Hubbard system suitable for topolog- 
ical quantum computing where fault-tolerance is a 
natural consequence. Alternatively, simulation of quan- 
tum systems by regime replication (where the Hamilto- 
nian of the simulator system is transformed to be equiv- 
alent to the Hamiltonian of some unknown system) is an 
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attractive application. 

Similar though modified dynamics found in other 
Hubbard-like systems should also be realizable in the 
atom-dressed photon superlattice, such as a glass phase 
stripes and other symmetry breaking phases [ssf . 
super-solid behavior [s^, etc. More importantly, since 
the excitation spectrum of a fluctuation is gapped in the 
Mott-insulator region, fluctuations due to temperature 
are exponentially suppressed. Thus, one will be able 
to observe a phase that is formally not a Mott insula- 
tor, but experimentally has very similar features lllllla . 
Disorder-induced destruction of the Mott state is 
also suppressed by this excitation-hole gap, and system 
tunability as via the Stark shift may allow for manual 
correction of these irregularities. Because of the Mott 
phase's robustness, devices based on this effect at non- 
zero temperature should be possible in this system. We 
also note that signatures of quantum many-body phe- 
nomena should appear in ID or finite arrays of cavities. 

Note Added: Whilst we were preparing our 
manuscript, we became aware of two related works, one 
by Hartmann et al. 37] , and one by Angelakis et al. |38l | , 
treating arrays of coupled quantum cavities. The former 
|87l | considers a four-state scheme for realising the re- 
quired quantum nonlinearity with time domain analysis. 
The latter [s^ l performs a similar analysis explicitly con- 
sidering two-state systems per cavity in a 1-D arrange- 
ment to realise an XY chain. Neither paper demonstrates 
Mott lobes as done here, nor did they perform a mean 
field analysis. As such, both works are qualitatively dif- 
ferent from the present one, but demonstrate that pho- 
tonic systems are ripe for the exploration of condensed 
matter physical effects. 



METHODS 



The eigenvalues and eigenvectors of the Hamiltonian 
of Eq.Elare 

^ |.g,n)+/3V^|e,n-l) ^ ^ 

(6) 



with eigenvalues 



S|±,„) = nio±x{n)-A/2. 



(7) 



Because n > 1, the above definition for the dressed states 
does not extend to n = 0, and so we must define the 
ground state for the dressed state system as |5,0), with 
E\a,o) = 0. 

To calculate the ground state wavefunction of our ex- 
tended Jaynes-Cummings Hamiltonian we must calculate 
the ground-state energy Eg as a function of -0 as follows. 
First: obtain the matrix elements of H^^^ in the onsite 
basis {\g, 0), |e, 0), \g, 0), |e, 1), . . .} truncated at a finite 
value of n states, Umax- Second: diagonalize this matrix 
(in the truncated basis) and identify the lowest eigen- 
value. Eg. Third: until Eg converges to 
its value at rimax = oo. Finally: minimize Eg with re- 
spect to V' for different values of k, lo, and A to obtain the 
phase diagram. At T = this approach is operationally 
equivalent I?! to the Gutzwiller Ansatz (GA) variational 
wave function technique often used, but can be 
extended to finite temperatures. 



For clarity, we will write out the first few terms of the mean-field Hamiltonian in the onsite basis, jg, 0), |e, 0), jg, 1), 



etc. 
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(8) 



where the entries in the boxes correspond to the usual block diagonal form of the Jaynes-Cummings Hamiltonian with 
the addition of the superfluid order parameter, and the number of excitations in each block increases by one as we 
move diagonally from block to block from the ground state \g, 0). The critical chemical potentials can be immediately 
determined by equating the lower eigenstates (negative dressed-state branch) from neighbouring blocks, with k = 0, 
which gives the result in Eq. [3 
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FIG. 1: (a) Schematic showing a two-dimensional array of photonic bandgap cavities, with each cavity containing a single 
two- level atom (spheres). The photonic bandgap is described by a periodic perturbations of the dielectric medium, in this 
case represented by the regions drilled through a thin membrane, for example via focussed ion beam. The spacing between 
the holes will be of order \/rj where A is the optical wavelength and rj the refractive index. The membrane should be of order 
\/{2r]) thick to confine light in the plane. Cavities are defined by defects in the periodic modulation, i.e. where a rod has not 
been drilled. By introducing a periodic array of defects in the photonic bandgap structure we can realise a lattice of photonic 
bandgap cavities, i.e. a photonic bandgap cavity superlattice. Photons can hop from one cavity to either of the three nearest 
neighbours with a frequency (photon hopping rate) k, illustrated for some of the cavities by the double headed arrows. To 
effect photon-photon repulsion, we insert a single two-state system into each cavity, represented here by the red spheres. If the 
membrane were fabricated in diamond, the two-state systems could be ion implanted nitrogen-vacancy (NV) centres. In this 
lattice we have chosen for each cavity to have three nearest neighbours, so that the regime where photonic repulsion dominates 
over hopping is easier to achieve, (b) Schematic showing some of the pertinent energy scales within each atom-cavity system. 
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FIG. 2: Eigenspectrum for a single atom in a high Q cavity, as a function of the atom-cavity detuning, centered around 
0, i.e. we plot E — nu. The eigenspectrum splits naturally into two branches, corresponding to the dressed states, \~\-,n) 
(upper branch) and | — , n) (lower branch). Note that the ground state \g, 0), naturally appears to be a member of each branch, 
introducing a significant departure from usual Bose-Hubbard style dynamics within each branch. Excluding 0), the branches 
anti-cross at A = 0, with the splitting increasing with increasing excitation number, n, given by the usual Rabi frequency, 
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Detuning, A [P] 

FIG. 3: Boundaries between Mott lobes in the limit of low tunneling (small k) as a function of jj, and A. The two lowest 
stable regions correspond to \g,0) on the left, and | — , 1) on the right. The boundary between these two domains follows the 
line jj, = -E|_j). Higher order states can be seen in the central region, defined by the critical chemical potentials fj,c{n) and 
correspond to | — ,2), | — ,3) etc. 
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FIG. 4: Slices showing the superfluid order parameter, i/), as a function of the photon hopping frequency, k, and the chemical 
potential, fi, for (a) detuning, A = 0, (b) A — —2, and (c) A = 2. The diagrams show Mott insulator lobes, indicated by the 
regions of i/; = 0, where the lowest have been shown. Dominating the left hand edge (where photonic repulsion dominates over 
hopping) is the Mott insulator phase (denoted MI), and the superfluid phase is found on the right hand edge (denoted SF). 
The white contour in (a) corresponds to the region where becomes nonzero, delineating the quantum phase transition. 
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FIG. 5: Plateaus with constant density, p, indicating regions with a definite number state excitations, as a function of the 
chemical potential, /i, and photon hopping frequency, k, for detuning, A = 0. The first three plateaus are indicated, and the 
ground state configurations correspond to |p, 0), | — , 1) and | — , 2) for 0, 1, and 2 excitations (photons in this case) respectively, 
with the plateaus shrinking in size with increasing excitation number. Regions with varying p have coherent states as the 
ground state configuration. 



